Observation of photon-assisted tunneling in optical lattices 
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We have observed tunneling suppression and photon-assisted tunneling of Bose- Einstein conden- 
sates in an optical lattice subjected to a constant force plus a sinusoidal shaking. For a sufficiently 
large constant force, the ground energy levels of the lattice are shifted out of resonance and tunnel- 
ing is suppressed; when the shaking is switched on, the levels are coupled by low-frequency photons 
and tunneling resumes. Our results agree well with theoretical predictions and demonstrate the 
usefulness of optical lattices for studying solid-state phenomena. 
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A number of experiments in recent years have shown 
that Bose-Einstein condensates (BECs) 1] loaded into 
optical lattices are well suited to simulating solid state 
systems 0, [3- Optical lattices are created by cross- 
ing two or more laser beams, and the resulting periodic 
potential landscapes (arising from the ac-Stark shift ex- 
erted on the condensate atoms) are intrinsically defect- 
free, their lattice wells have controllable depths, and it is 
possible to move or accelerate the entire structure. This 
flexibility has made it possible to study dynamical effects 
such as Bloch oscillations Q and resonant tunneling Q 
as well as ground-state quantum properties such as the 
Mott-insulator transition [gi] . More recently, the coherent 
suppression of inter-well tunneling by strong driving of 
the lattice has been demonstrated In this Letter, we 
explore an effect Q that is analogous to photon-assisted 
tunneling in solids and arises from the interplay between 
static acceleration and strong driving of the lattice. We 
observe two regimes, a linear and a nonlinear one, with 
different dependencies of the observed tunneling on the 
theoretically predicted behaviour. 

Photon-assisted tunneling occurs when adjacent poten- 
tial wells whose ground states are tuned out of resonance 
by a static potential are coupled by photons (see Fig. 1). 
The static force leads to a suppression of resonant tunnel- 
ing between the ground states. This suppression and the 
related Wannier-Stark localization of the wavefunction 
have been intensively discussed in the theoretical litera- 
ture In this work, we report a direct measure- 
ment of this suppression based on the spatial tunneling 
of the condensate atoms. When photons of an appro- 
priate frequency are present whose energy bridges the 
gap created by the static potential, tunneling is (partly) 
restored. In solid state systems, the photons are typ- 
ically in the microwave frequency range and the static 
potential is provided by an electric bias field applied to 
the structure. So far, photon-assisted tunneling has been 
observed in superconducting diodes semiconductor 
superlattices [ij, [H, [l^ a^nd quantum dots [l^, [lB| ■ 

Our system consists of a BEG inside a one-dimensional 
optical lattice. The static potential is provided by a con- 
stant acceleration of the lattice, resulting in a constant 
force F in the lattice rest frame and hence in a potential 
difference AE' = Fdi, between adjacent wells a distance 




FIG. 1: Photon-assisted tunneling, (a) In an optical lattice at 
rest, the ground-state levels are resonantly coupled, leading to 
a tunneling energy J. (b) When a linear potential is applied, 
e.g. by accelerating the lattice, the levels are shifted out of 
resonance and tunneling is suppressed, (c) If the lattice is now 
periodically shaken at an appropriate frequency, the levels can 
again be coupled through photons of energy fiu) and tunneling 
is partially restored. 



apart. The role of the photons is played by a periodic 
shaking of the lattice at frequency oj that leads to the 
creation of sidebands around the carrier frequency of the 
laser beam. In the limit of sufficiently deep lattice wells 
and neglecting higher-lying energy levels, our system can 
be described by the Hamiltonian [Sj 

+ /S.E jfij + K cos(cjt) jfij , 
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where (fp are the boson creation and annihilation oper- 
ators on site z, hi = c\ci are the number operators, and 
K and oj are the strength and angular frequency of the 
shaking, respectively. The first line of this equation is 
the Bose-Hubbard model 17[ with the tunneling matrix 
element J and the on-site interaction term U (in a BEC, 
the on-site interaction is due to atom-atom collisions and 
hence proportional to the s-wave scattering length and 
the density of the BEC). In the second line, the first term 
describes the constant potential, whereas the second term 
represents the sinusoidal shaking of the lattice. While for 
a sufficiently strong linear potential inter-well tunneling 
is suppressed, leading to Wannier-Stark localization, re- 
cent theoretical work [8] predicts that the shaking term 
can partially restore it, leading to an effective tunneling 
rate 



\J,s{K„)/J\ = \Jn{Ka)\ 
when the resonance condition 

nhiij — FcIl 



(2) 



(3) 



is satisfied, where n is an integer denoting the order of the 
photon-assisted resonance, J7n is the n-th order ordinary 
Bessel function, and Kq = K/fuv is the dimensionless 
parameter characterizing the shaking amplitude. 

In our experiment we produced BECs of ^^Rb con- 
taining around 5 x 10* atoms in a crossed optical dipole 
trap. The two dipole traps were created by gaussian laser 
beams at 1030 nm wavelength and a power of around 1 W 
per beam focused to waists of 50 fim, and the frequen- 
cies of the resulting trapping potentials could be con- 
trolled independently. Subsequently, the BECs held in 
the dipole trap were loaded into an optical lattice cre- 
ated by two counter-propagating gaussian laser beams 
(A = 852 nm) with 120 /xm waists by ramping up the 
power of the lattice beams in about 50 ms. The resulting 
periodic potential V{x) = Vq sin'^(7rx/c?i,) had a lattice 
spacing — A/2 = 426 nm and its depth Vq was mea- 
sured in units of the recoil energy E^^c = tt^ / {2md\) , 
where m is the mass of the Rb atoms. By introducing 
a frequency difference Ai/ between the two lattice beams 
(using acousto-optic modulators, which also control the 
power of the beams) , the optical lattice could be moved 
at a velocity v = d^Av or accelerated with an accel- 
eration a = dL^^^- In order to periodically shake the 
lattice, Av was sinusoidally varied with frequency uj and 
amplitude At'max leading to a time-varying force (in the 
rest frame of the lattice) 

F{t) = mudLAvra^-^ coslijjt) = Finax cos{ujt) . (4) 

The dimensionless shaking parameter is then given 
by 

= K/huj = md\Al'^^y,/fl = TT^ Al/max/2t^rcc- (5) 

Our method for measuring the effective tunneling pa- 
rameter I Jeff/"/ 1 (where the modulus indicates that this 



measurement is not sensitive to the si gn o f J^s) is based 
on the free expansion of the BEC 0,13 confined only 
radially but free to move along the direction of the lattice 
(the dipole trap frequency in that direction being on the 
order of a few Hz and hence negligible for our purposes). 
After condensation was reached in the crossed dipole trap 
and the optical lattice had been ramped up, the trapping 
beam perpendicular to the lattice direction was suddenly 
switched off. Subsequently, the in-situ width of the BEC 
in the lattice direction was measured by flashing on a 
resonant beam and imaging the shadow cast by the BEC 
on a CCD camera. In the experiments with an acceler- 
ated and / or shaken lattice, | Jeff/J| was determined by 
measuring the expansion for the same lattice depth both 
in the driven case and without the driving. 

In a preliminary experiment, we studied the tunnel- 
ing suppression caused by shifting adjacent ground states 
out of resonance through a constant force F acting on 
the BEC inside a lattice that was subjected to an ac- 
celeration a. Typical accelerations in our experiment 
were between and 4ms~^, meaning that for an expan- 
sion time of 400 ms the lattice (and therefore the BEC) 
would have been displaced by up to 32 cm, two orders 
of magnitude more than the field-of-view of our imag- 
ing system. Also, a displacement of more than 100 /im 
along the lattice direction would have led to a restoring 
force of the longitudinal harmonic trap created by the 
radial dipole trap beam corresponding to an acceleration 
arcstorc > 0.1 ms~^. Therefore, in order to be able to 
achieve a high resolution in our measurements of the ex- 
pansion rate (and hence J) of the condensate, implying 
a long expansion time, whilst keeping the displacement 
of the lattice below « 100 /im, we used a rectangular 
acceleration profile that alternated between +a and —a 
and therefore 'rocked' the lattice back and forth. In this 
way, the modulus of the resulting force and hence the 
energy shift between adjacent wells was constant, while 
the mean position of the lattice (and hence the BEC) 
remained close to the center of the dipole trap. In or- 
der to separate the frequency regimes of the rocking and 
the shaking motion, we chose a 'rocking frequency' of 
30 Hz, which was much smaller than the resonant fre- 
quencies for photon-assisted tunneling (typically around 
150 - 400 Hz). 

Figure 2 (a) shows the results of our measurements 
of I Jeff I in an accelerated (rocked) lattice. As expected, 
when the energy difference FdL between adjacent lev- 
els is increased, resonant tunneling is reduced and, for 
FdL ~ J, completely suppressed (as recently also ob- 
served for single-atom tunneling in a double-well struc- 
ture [l^). In this limit, the energy levels in the individual 
wells can be viewed as Wannier-Stark levels. Our data 
are fitted very well by a Lorentzian, but to our knowl- 
edge there is no analytical prediction for such a depen- 
dence in the theoretical literature. In [I^Ij an expression 
for the wavefunction in a tilted lattice as a function of 
the applied force is given, but an analytical calculation 
of Jcs{F) has yet to be done. 
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FIG. 2: (a) Suppression of tunneling by a linear potential. 
Shown here is the normalized effective tunneling parameter 
\JcS / J\ as a function of the linear potential FdL in units of 
the tunneling energy. When Fdhj J « 1, the ground state lev- 
els are shifted out of resonance and tunneling is suppressed 
almost completely. The solid line is a Lorentzian fit with a 
half-maximum half- width of 0.13. (b) Photon-assisted tunnel- 
ing resonances in a shaken lattice. For a fixed linear potential 
Fdhj h — 380 Hz, the condensate width after 400 ms of free ex- 
pansion is plotted as a function of the normalized shaking fre- 
quency hw/Fdh- The fixed shaking parameter was Ko = 1.8 
for the one-photon resonance (solid circles) and Kq =3.1 for 
the two-photon resonance (open circles), corresponding to the 
first maximum of the and jTij Bessel functions, respectively. 
For both graphs, Vo/E^oc = 5 and J/h = 380 Hz. 



Tunneling between the on-site levels shifted out of res- 
onance by the static acceleration can be partially restored 
by sinusoidally shaking the lattice at a frequency uj sat- 
isfying the resonance condition of Eq. [3] The shaking 
of the lattice effectively creates low-frequency 'photons' 
that bridge the energy gap between adjacent wells with 
n such photons (see Fig. 1 (c)). Figure 2 (b) shows 
the condensate width after 400 ms of free expansion in- 
side a rocked lattice with FcIl/J = 1 as a function of 
the shaking frequency uj. One clearly sees two photon- 
assisted tunneling peaks at huj = FcLl (for Kq = 1.8, 
where Si{Kq) has its first maximum) and at 2fiLo = Fd^ 
(for Kq = 3.1, the first maximum of J2{Kq)). The peaks 



are extremely narrow, with a width of « 3 Hz (when we 
repeated the experiment with smaller values of FcLl, this 
width increased sHghtly). The two side-peaks of each 
peak are evidence of additional photon- assisted tunnel- 
ing events due to the rocking motion of the lattice (for 
the two-photon resonance at 2hLjj = FcIl, the side-peaks 
are at half the distance to the main peak compared to 
the one-photon expected) . 
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FIG. 3: Photon-assisted tunneling as a function of the shaking 
parameter Ka. Shown here are the one-photon resonance at 
ijj/2-K — 380 Hz (a) and the two-photon resonance at uj/2n = 
190 Hz (b). In both graphs, the full and open squares are the 
measurements for A*' ~ 5 x lO** and A*' ~ 0.5 x lO** , respectively. 
The solid lines are the moduli of the Ji{Ko) and J2{Kq) 
Bessel functions, respectively, whereas the dashed lines are 
the squares of these functions. The lattice depth Vo/E-rcc = 5 
and the constant force FdL/h — 380 Hz in this experiment, 
with a free expansion time t — 100 ms. 

Finally, we studied the dependence of the effective 
tunneling rate on the shaking parameter Kq- Figure 
3 summarizes our results for the one-photon and two- 
photon resonances. Theory predicts [8] that the effec- 
tive tunneling rates \Jcft/J\ for resonances of rt-th order 
should vary as an n-th order ordinary Bessel function (see 
Eq. [2) . While qualitative agreement between experiment 
and theory is good, with the positions of the maxima 
and minima of the one- and two-photon resonances as a 
function of Kq coinciding perfectly with theoretical pre- 
dictions, the absolute values of \Jcs/J\ he consistently 
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below the theoretical curves by a factor of about 1.3. 

Interestingly, quantitative agreement between experi- 
ment and theory is better if we use the squares of the 
Bessel functions rather than their moduli. A dependence 
of the photon-assisted tunneling rate on the square of 
the Bessel function is expected, e.g., for experiments on 
Josephson junctions irradiated by microwaves and, more 
generally, if the tunneling in a multi-well structure is se- 
quential rather than coherent. The difference between 
coherent and sequential tunneling has been extensively 
studied in the theoretical literature [131 • Sequential tun- 
neling would require a dephasing mechanism between two 
successive tunneling events. In our experiments, such a 
mechanism could be the dynamical instability inside the 
optical lattice [U [l^, [23]: owing to the constant ap- 
plied acceleration a, the EEC moves through the Bril- 
louin zone and hence through the dynamically unstable 
region within the Brillouin zone between about 0.5proc 
(where Proc = ttH/cIl is the recoil momentum ) and the 
zone edge. For the values of the intrinsic nonlinearity of 
the BEC in our experiments (determined by the number 
of atoms and the trap frequencies), we found in previous 
experiments (2l| (and have verified for this work) that 
the BEC loses its phase coherence after a few millisec- 
onds inside the unstable region, which is comparable to 
the Bloch periods in the present experiment. Since the 
effective tunneling frequencies JoS in our experiment are 
less than 100 Hz, the corresponding dephasing rate is al- 
most an order of magnitude larger and hence it is likely 
that dephasing of neighbouring wells occurs between two 
tunneling events. 

In future experiments one might use, e.g., Feshbach 
resonances in order to tune the nonlinearity and hence 
move from the strongly nonlinear regime to the linear 
regime in order to test our hypothesis. As a prelimi- 
nary test, we have repeated the measurement of the one- 
photon assisted tunneling rate as a function of Kq with 
the smallest atom number that allowed us to measure the 
free expansion {N « 0.5 x 10"*), resulting in a condensate 
density that was about a factor of 3 smaller. Again, 
we obtained qualitatively similar results to the measure- 
ments with « 5 X 10'*, but this time the absolute values 
of \Jcff/J\ agreed better with the linear Bessel- function 



prediction (see Fig. 3). 

A further indication that dephasing might be respon- 
sible for the observed deviation from the linear Bessel 
functions comes from our measurements of the dynami- 
cal suppression of tunneling in a shaken lattice without 
linear acceleration Q, in which our measured values for 
Jeff/ J agreed perfectly with the linear Bessel function 
prediction. In that system, the BEC does not cross the 
unstable region of the Brillouin zone, and we have ex- 
perimentally verified that during the shaking the BEC 
retains its phase coherence. 

If the two regimes of photon-assisted tunneling ob- 
served in our experiment do, indeed, correspond to coher- 
ent and sequential tunneling, our system is ideally suited 
to studying the cross-over between these two extremes 
in a well-controlled way. However, we also have to con- 
sider the possibility that other effects play a role. In 
particular, it is conceivable that a self-trapping mecha- 
nism is present that depends on the relative magnitude 
of the nonlinearity U and the effective tunneling param- 
eter Jeff. Self-trapping in static optical lattices has al- 
ready been observed [Isl] and could, if present also in our 
strongly driven system, lead to a suppression of the tun- 
neling rate compatible with our observations, also for the 
case of a static force as in Fig. 2(a). On the other hand, 
our system is only in a self-trapping regime (as calculated 
from U/h w 10 — 30 Hz and Jeff) in a small region around 
the zeroes of the Bessel functions, whereas we observe the 
largest deviation from the linear prediction close to the 
local maxima. 

In summary, we have demonstrated photon-assisted 
tunneling of BECs in a linearly accelerated and sinu- 
soidally shaken optical lattice. Our results agree quan- 
titatively with recent theoretical predictions and show 
the need for a theoretical investigation into the difference 
between the (roughly) linear and the nonlinear regimes 
for which we observe difference dependencies of the effec- 
tive tunneling rate on the predicted linear Besse-function 
scaling. 
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